Chapter 1 

Positive energy theorems in 
General Relativity 

1.1 Introduction 

The aim of this chapter is to present an introduction and also an overview 
of some of the most relevant results concerning positivity energy theorems in 
General Relativity. These theorems provide the answer to a long standing prob- 
lem that has been proved remarkably difficult to solve. They constitute one of 
the major results in classical General Relativity and they uncover a deep self- 
consistence of the theory. 

In this introductory section I would like to present the theorems in a complete 
form but with the least possible amount of technical details, in such a way that 
the reader can have a rough idea of the basic ingredients. The examples that 
illustrate the hypothesis of the theorems are discussed in the following sections. 

An isolated system is an idealization in physics that assumes that the sources 
are confined to a finite region and the fields are weak far away from the sources. 
This kind of systems are expected to have finite total energy. In General Rela- 
tivity there are several ways of defining isolated systems. For our purpose the 
most appropriate definition is through initial conditions for Einstein equations. 
The reasons for that are twofold. First, the notion of total energy has been 
discovered and formulated using a Hamiltonian formulation of the theory which 
involves the study of initial conditions. We refer the reader to the chapter of 
Domenico Giulini for this topic. Second, the proofs of the positive mass the- 
orem are mainly given in terms of initial conditions. For a discussion of the 
initial value formulation of Einstein equations we refer to the chapter of James 
Isenberg. 

Initial conditions for Einstein equations are characterized by initial data set 
given by (S, hij, Kij, /i, j l ) where S is a connected 3-dimensional manifold, a 
(positive definite) Riemannian metric, K%j a symmetric tensor field, \i a scalar 
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field and f a vector field on S, such that the constraint equations 

DjK ij - D i K = -8irf, (1.1) 
R - K l0 K l] +K 2 = I6717/, (1.2) 

are satisfied on S. Here D and R are the Levi-Civita connection and scalar 
curvature associated with hij, and K = Kijh lJ . In these equations the indices 
i,k, . . . are 3-dimensional indices, they are raised and lowered with the metric 
hij and its inverse h lJ . The matter fields are assumed to satisfy the dominant 
energy condition 

/' \ T 7 ^- (1.3) 

The initial data model an isolated system if the fields are weak far away from 
sources. This physical idea is captured in the following definition of asymptoti- 
cally flat initial data set. Let Br be a ball of finite radius R in R 3 . The exterior 
region U — M 3 \ Br is called an end. On U we consider Cartesian coordinates 

( 3 ■ \ 1 / 2 

x l with their associated euclidean radius r — [^2i—i( x Jj an( i ^ &ij b° 
the euclidean metric components with respect to i 1 . A 3-dimensional manifold 
S is called Euclidean at infinity, if there exists a compact subset K, of S such 
that S \ JC is the disjoint union of a finite number of ends U k - The initial data 
set (S,hij, Ktj, /i, j l ) is called asymptotically flat if S is Euclidean at infinity 
and at every end the metric hij and the tensor Kij satisfy the following fall off 
conditions 

= 6ij + 7y, = 0{r- 2 ), (1.4) 

where 7y = O^ 1 ), d k jij = 0(r" 2 ), did k Jij = 0(r" 3 ) and d k Kij = 0(r" 3 ). 
These conditions are written in terms of Cartesian coordinates x l attached at 
every end Uk- Here di denotes partial derivatives with respect to these coordi- 
nates. 

At first sight it could appear that the notion of asymptotically flat manifold 
with "multiple ends" Uk is a bit artificial. Certainly, the most important case 
is when S = K 3 , for which this definition trivializes with K, = Br and only 
one end U = M 3 \ Br. Initial data for standard configurations of matter like 
stars or galaxies are modeled with S = M 3 . Also, gravitational collapse can 
be described with this kind of data. However, initial conditions with multiple 
ends and non-trivial interior K, appear naturally in black hole initial data as we 
will see. In particular, the initial data for the Schwarzschild black hole has two 
asymptotic ends. On the other hand, this generalization does not imply any 
essential difficulty in the proofs of the theorems. 

Only conditions on hij and are imposed in (|1.4p and not on the matter 
fields /j, and j l , however since they are coupled by the constraint equations 
1|1.1|) - 1|1.2|) the fall off conditions (|1.4j) impose also fall off conditions on 

The fall off conditions (|1.4p are far from being the minimal requirements for 
the validity of the theorem. This is a rather delicate issue that have important 
consequences in the definition of the energy. We will discuss this point in section 
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11.21 We have chosen these particular fall off conditions because they are simple 
to present and they encompass a rich family of physical models. 

For asymptotically flat initial data the total energy and linear momentum 
of the spacetime are defined as integrals over 2-spheres at infinity at every end 
by the following expressions 



where s l is its exterior unit normal and dso is the surface element of the 
2-sphere with respect to the euclidean metric. We emphasize that for every end 
Uk we have a corresponding energy and linear momentum E^,PJ k -., which can 
have different values. We will discuss examples of that in section IL21 

The quantities E and Pi are defined on the asymptotic ends and they depend 
only on the asymptotic behaviour of the fields and Kij. However, since htj 
and satisfy the constraint equations (|1.1|) - (|1.2[) and the dominant energy 
condition (ll.3[) holds these quantities carry in fact information of the whole 
initial conditions. 

The energy E and the linear momentum Pi are components of a 4-vector 
P a = (E,Pi) (indices a,b, c, ... are 4-dimensional). We will discuss this in 
section 11.31 The total mass of the spacetime is defined by 



We have all the ingredients to present the positive energy theorem. 

Theorem 1.1.1 (Positive ene rgy theorem). Let (S, hij, Kij, pi,j l ) be an asymp- 
totically flat (with possible many asymptotic ends), complete, initial data set, 
such that the dominant energy condition hl.3\) holds. Then the energy and lin- 
ear momentum (E,Pi) defined by i U.5)) - (|1.6p satisfies 



at every end. Moreover, E = at any end if and only if the initial data corre- 
spond to the Minkowski space-time. 

The word "complete" means that (S,hij) as Riemannian manifold is com- 
plete. That is, no singularities are present on the initial conditions. But the 
space-time can be singular since singularities can developed from regular initial 
conditions, for example in the gravitational collapse. We will discuss that in 
more detail in section [i~2"1 

Note that Theorem 11.1. II allows the vector P a to be null and non trivial. 
However, it has been shown in [2] that if the energy momentum vector P a is 
null then it vanishes identically. 




(1.5) 



(1.6) 




(1.7) 




(1.8) 
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One remarkable aspect of this theorem is that it is non-trivial even in the case 
where S = M. 3 and no matter fields /i = j l = are present. This correspond to 
the positivity of the energy of the pure vacuum gravitational waves. We present 
explicit examples of this in section ITT21 

For spacetimes with black holes there are spacelike surfaces that touch the 
singularity. For that kind of initial conditions theorem II. 1.11 does not apply. 
Physically it is expected that it should be possible to prove a positivity energy 
theorem for black holes without assuming anything about what happens inside 
the black hole. That is, it should be possible to prove an extension of the positive 
energy theorem for initial conditions with inner boundaries if the boundary 
represents a black hole horizon. The following theorem deals precisely with 
that problem. 

Theorem 1.1.2 (Positive energy theorem with black hole inner boundaries). 

Let (S,hij,Kij) be an asymptotically flat, complete, initial data set, with S = 
M 3 \ B, where B is a ball. Assume that the dominant energy condition (|1.3|) 
holds and and that dB is a black hole boundary. Then the energy momentum 
E, P l defined by \1.5\) - (|1.6|) satisfies 

E > \/P l P % > 0. (1.9) 

Moreover, E = if and only if the initial data correspond to the Minkowski 
space-time. 

We will explain what are black hole inner boundary conditions in section 

PI 

The plan of the chapter is the following. In section[0]we discuss the concept 
of the energy E and we present examples that illustrate the hypothesis of the 
positive energy theorem. In section [TT51 we analyze the linear momentum Pj and 
describe its transformation properties. In section [L4l we review the main steps of 
the proof of theorems II . 1.11 and 11.1.21 Finally in section [T3] other recent related 
results are discussed and the relevant current open problems are presented. 

1.2 Energy 

A remarkable feature of the asymptotic conditions (|1.5p is that they imply 
that the total energy can be expressed exclusively in terms of the Riemannian 
metric hij of the initial data (and the linear momentum in terms of hij and the 
second fundamental form Kij). Hence the notion of energy can be discussed 
in a pure Riemannian setting, without mention the second fundamental form. 
Moreover, as we will see, there is a natural corollary of the positive energy 
theorem for Riemannian manifolds. This corollary is relevant for several reasons. 
First, it provides a simpler and relevant setting to prove the positive energy 
theorem. Second, and more important, it has surprising applications in other 
areas of mathematics. Finally, to deal first with the Riemannian metric and 
then, in the next section, with the second fundamental form to incorporate 
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the linear momentum, reveal the different mathematical structures behind the 
energy concept. 

In the previous section we have introduced the notion of an end U, the energy 
is defined in terms of Riemannian metrics on U. To emphasize this important 
point we isolate the notion of energy defined in the introduction in the following 
definition. 

Definition 1.2.1 (Energy). Let hij be a Riemannian metric on an end U given 
in the coordinate system x % associated with U . The energy is defined by 



Note that in this definition there is no mention to the constraint equations 
(|l.l[) - (|1.2p . Also, the definition only involve an end U, there is no assumptions 
on the interior of the manifold. 

In the literature it is custom to call E the total mass and denote it by m or 
M . In this article, in order to emphasize that E is in fact the zero component 
of a four vector we prefer to call it energy and reserve the name mass to the 
quantity M defined by (|1.7p . When the linear momentum is zero, both quantities 
coincides. 

The definition of the total energy has three main ingredients: the end U, 
the coordinate system x l and the Riemannian metric hij. The metric is al- 
ways assumed to be smooth on U, we will deal with singular metrics but these 
singularities will be in the interior region of the manifold and not on U. 

There exists two potentials problems with the definition 11.2.11 The first 
one is that the integral (|1.10j) could be infinite. The second, and more subtle, 
problem is that the mass seems to depend on the particular coordinate system 
x 1 . Both problems are related with fall off conditions for the metric. In the 
previous section we have introduced in equation (|1.4[) an example of this kind 
of conditions. These conditions are probably sufficient to model most physically 
relevant initial data. However, it is interesting to study the optimal fall off 
conditions that are necessary to have a well defined notion of energy and such 
that the energy is independent of the coordinate system. 

To study this problem, we introduce first a general class of fall off conditions 
as follows. Given an end U with coordinates x % , and an arbitrary real number 
a, we say that the metric hij on U is asymptotically flat of degree a if the 
components of the metric with respect to these coordinates have the following 
fall off in U as r — > oo 



with 7y = 0(r~ a ), dkjij = 0(r _Q_1 ). The subtle point is to determine the 
appropriate a decay. To understand the meaning of this coefficient let us discuss 
the following relevant example given in [19] (see also [6]). Take the euclidean 
metric Sij in Cartesian coordinates x % and consider coordinates y % defined by 




(1.10) 



(l.n) 



(1.12) 
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where p is defined by 



r = 



p + cp 



l-a 



(1.13) 



for some constants c and a. Note that p = ^X)i=i(y*) 2 j ■ The components 
g[ ■ of the euclidean metric in coordinates y % have the following form 



9. 



(1.14) 



where 7^ satisfies the decay conditions (|l.lip with the arbitrary a prescribed 
in the coordinate definition (|1.13[) . That is, the metric in the new coordinate 
system y l is asymptotically flat of degree a. 

We calculate the energy in the coordinates y l using the definition (|1.10|) . We 
obtain 



Of course, we expect that the energy of the euclidean metric should be zero in 
any coordinate system. The interesting point of this example is the limit case 
a = 1/2, the example shows that if the energy has any chance to be coordinate 
independent, then we should impose a > 1/2. The following theorem, proved 
in [3] and [12], says that this condition is also sufficient. 

Theorem 1.2.2. Let U be an end with a Riemannian metric hij such that is 
satisfies the fall off conditions (|1.11|) with a > 1/2. Assume also that the scalar 
curvature R is integrable in U , that is 



Then the energy defined by (jl.lOp is unique and it is finite. 

In this theorem unique means if we calculate the energy in any coordinate 
system for which the metric satisfies the decay conditions (jl.lip with a > 1/2 
we obtain the same result. This theorem ensure that the energy is a geometrical 
invariant of the Riemannian metric in the end U. Historically, this theorem was 
proved after the positive energy theorems. In the original proofs of the positive 
energy theorems different decay conditions for the metric have been used. The 
decay conditions are usually formulated in terms of integrals of derivatives (i.e. 
Sobolev spaces) (see [3]) which are more flexible for many applications. This 
particular formulation (which is simpler to present) of theorem 11.2.21 was taken 
from [T3]. The decay conditions with a > 1/2 together with the condition p. 161) 
on the scalar curvature are called mass decay conditions. The freedom in the 
coordinates x % is only a rigid motion at infinity (see [3]). 

Theorem 11.2 . 21 completes the geometric characterization of the energy at the 
end U. We turn now to positivity. It is clear that the energy can have any sign 




00, a < 1/2, 
c 2 /8, a =1/2, 
0, a>l/2. 



(1.15) 




(1.16) 
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on U. The model example is given by the initial data for the Schwarzschild 
black hole, with metric on U given by 

hij = il>%, (1.17) 

where ip is the following function 

^ = l + g, (1.18) 

with C an arbitrary constant. Computing the energy for this metric we obtain 
E = C. The constant C can of course have any sign. It is however important 
to emphasize that theorem 11.2.21 asserts that the energy is well defined and it is 
an invariant of the geometry of the end even when it is negative. 

To ensure the positivity of the energy we need to impose two important 
conditions. One is a local condition: the positivity of the local energy given 
by the dominant energy condition (ll.3|) . The other is a global condition on the 
manifold: the manifold should be complete or should have black hole boundaries. 

Initial conditions with 

Kij = 0, (1.19) 

are called time symmetric initial data. That is, time symmetric initial data 
are characterized only by a Riemannian metric hij. Conversely, any Rieman- 
nian metric can be interpreted as a time symmetric initial data. However, an 
arbitrary metric will not satisfy the dominant energy condition (|1.3j) . In ef- 
fect, inserting condition (|1.19j) in the constraint equation (|1.2|) and using the 
dominant energy condition (|1.18[) we obtain 

R > 0. (1.20) 

Only metrics that satisfy (| 1 . 20[) can be interpreted as time symmetric initial 
data for which the dominant energy condition holds. But then, any metric 
such that (ll.20[) holds satisfies the dominant energy condition and it is a good 
candidate for the positive energy theorem. And hence we obtain the following 
corollary of theorem II. 1.11 

Corollary 1.2.3 (Riemannian positive mass theorem). Let (S,hij) be a com- 
plete, asymptotically flat, Riemannian manifold. Assume that the scalar curva- 
ture is non-negative (i.e. condition (jl -20[) ). Then the energy is non-negative at 
every end and it is zero at one end if an only if the metric is flat. 

This corollary was proved with the optimal decay conditions for the metric 
in [3] and [2"6] . 

The interesting mathematical aspect of this corollary is that there is no 
mention to the constraint equations, the second fundamental form or the matter 
fields. This theorem is a result in pure Riemannian geometry. It has surprising 
applications in the solution of the Yamabe problem (see the review article [26] 
and reference therein). 
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Note that it is not necessary to impose that the whole second fundamental 
form is zero to have (ll.20p , from equation (|1.2p it is clear that is enough to have 
K = 0. This class of initial data are called maximal and they have important 
properties (see the chapter by J. Isenberg). In particular, positive energy theo- 
rems for this kind of data are easier to prove (mainly because condition (ll.20|) 
holds) than for general initial data. 

Let us discuss some examples of corollary 11.2.31 We begin with the case 
with one asymptotic end and trivial topology, namely S = K 3 . For arbitrary 
functions ip, metrics of the form (jl.171) are called conformally flat, they provide 
a very rich family of initial conditions which have many interesting applications 
(for example, initial data for black hole collisions, see the review article [15]). 
The scalar curvature for this class of metrics is given by 

R = -8i>- 5 Ai), (1.21) 

where A is the euclidean Laplacian. If ip satisfies the fall off conditions 

tp = l + u, u = 0(r- 1 ), d k u = 0(r- 2 ), (1.22) 

then the energy for this class of metric is given by 

E= — - lim I d r ^dsQ. (1.23) 

For i/j given by (11.181) we obtain R — 0, and then the metric satisfies the 
local condition f| 1 . 20[) for any choice of the constant C. However, this metric 
can not be extended to R 3 since the function ip is singular at r = and hence, 
as expected, corollary 1 1 . 2 . 31 docs not apply to this case. Let us try to prescribe 
a function with the same decay (and hence identical energy) but such that it is 
regular at r = 0. For example 

i/> = ! + —= ¥■ . (1.24) 

Using (11.23)) we obtain again that E — C. For any value of C the function tp 
is strictly positive and bounded on R 3 and hence the metric is smooth on R 3 . 
That is, it satisfies the completeness assumption in corollarv ll.2.31 Using (|1.2ip 
we compute the scalar curvature 

We have R > if and only if G > 0. Also, in this example the mass is zero if 
and only if the metric is flat. 

Other interesting examples can be constructed with conformally flat metrics 
as follows. Let ip be a solution of the Poisson equation 



(1.26) 
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that satisfies the decay conditions (|1.22[) , where fx is a non- negative function of 
compact support in R 3 . Solution of (|1.26p can be easily constructed using the 
Green function of the Laplacian. By equation (|1.2ip , the scalar curvature of 
the associated conformal metric (|1.17|) will be non-negative and the function ft, 
is related to the matter density \x by 

M =-|-=/^- 5 . (1-27) 

Note that we can not prescribe, in this example, exactly the matter density 
fi, we prescribe a conformal rescaling of \x. However, it is enough to control de 
support of fi. The support of /i represents the localization of the matter sources. 
Outside the matter sources the scalar curvature (for time symmetric data) is 
zero. 

For conformally flat metrics in R 3 there is a very simple proof of corollary 
11.2.31 We write equation (|1.21|) as 

^5 ) <ip6 



Integrating this equation in R 3 , using for the first term in the right-hand side 
the Gauss theorem, the condition tjj — > 1 as r — > oo and the expression (|1.23|) 
for the energy we finally obtain 



2ir 7 R 3 V 8 ^ 

where dv$ is the flat volume element. This formula proves that for metric of the 
form l|1.17l) we have E > if R > and E = if and only if = . This 
proof easily generalize for conformally flat maximal initial data. 

Asymptotically flat initial conditions in R 3 with no matter sources (i.e. fi = 
f = 0) represent pure gravitational waves. They are conceptually important 
because they describe the dynamic of pure vacuum, independent of any matter 
model. Note that in that case the energy condition (jl.3|) is trivially satisfied. 

In the previous examples the only solution with pure vacuum R = in R 3 
is the flat metric, because by equation (|1.2ip we obtain Atp = and the decay 
condition (|1.22[) implies tp = 1. In order to construct pure waves initial data we 
allow for more general kind of conformal metrics, let be given by 



- 2q (dp 2 + dz 2 )+p 2 dip 2 ], (1.30) 



where (p, z, tp) are cylindrical coordinates in R 3 and the functions q and a depend 
only on (/?, z). That is, the metric hij given by (|1.30|) is axially symmetric. 
The scalar curvature of the metric (| 1 . 30[) is given by 

_ l Re (^) = l A o + ±\da\ 2 - Ja 29 , (1.31) 
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where A, as before, is the 3-dimensional fiat Laplacian and A2 is the 2-dimensional 
Laplacian in cylindrical coordinates given by 

A 2 q = d 2 p q + d 2 z q. (1.32) 

If we impose R = 0, equation (|1.31[) reduce to 

A^-^A 2 q = 0, (1.33) 

where ip 4 = e a . To construct metrics of the form (|1.30j) that satishcs R = a 
function q is prescribed and then the linear equation (|1.33|) is solved for ip. The 
function q can not be arbitrary, it should satisfy a global condition (which is 
related with the Yamabe problem mentioned above), see [TT] for details. This 
kind of metric are called Brill waves, they have been used by D. Brill in one of 
the first proofs of the positive energy theorem [5] . Let us discuss this proof. 

In order to be smooth at the axis the metric (|1.30[) should satisfies q = 
at p = 0. For simplicity we also impose a strong fall off condition on q at 
infinity, namely q — 0(r -2 ), diq = 0(r~ 2 ). For a we impose a = 0(r~ r ) and 
di<r = 0(r~ 2 ). Using these decay assumptions is straightforward to check that 
the energy of the metric (|1.30l) is given by 

E= — - lim * d r ads . (1.34) 

By Gauss theorem, using that q = at the axis and the fall off condition of 
q at infinity we obtain that 

A 2 qdv = 0. (1.35) 

Integrating equation (|1.31j) in R 3 , using (|1.35|) and using the expression 
(11.34|) for the energy we obtain 

e= ^L G |9ct|2 + Rea2q ) dv °' (L36) 

That is, R > implies E > 0. In particular for vacuum R = 0, we have 

E = TT I \M 2 dv Q . (1.37) 

This positivity proof can be extended in many ways, in particular it has ap- 
plications for the inequality between energy and angular momentum discussed 
in section 11.51 (see the review article [T7] and the lectures notes [13] [2] , and 
reference therein). 

We turn now to manifolds with many asymptotic flat ends and interior K, 
with non-trivial topology defined in scction ri.il Let us first present some basic 
example of the definition of asymptotic euclidean manifold, without mention 
the metric. 
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Take out a point in R 3 , the manifold S = M 3 \ {0} is asymptotic Euclidean 
with two ends, which we denote by Uq and Ui. In effect, let B 2 and B\ be two 
balls centered at the origin with radius 2 and 1 respectively. Define JC be the 
annulus centered at the origin B 2 \ B± . Then S \ JC has two components Uq and 
Ui, where U = K 3 \ B 2 and Ui=B 1 \ {0}. The set U is clearly an end. The 
set U\ is also an end since the a ball minus a point is diffeomorphic to M 3 minus 
a ball. This can be explicitly seen using Cartesian coordinates centered at the 
origin x l , then the map given by the inversion 

l/^r-V, (1.38) 

provide the diffeomorphism between M 3 \ B\ and B\ \ {0}. 

In the same way M 3 minus a finite number N of points ik is an Euclidean 
manifold with N + 1 ends. For each ik take a small ball Bk of radius rrjA , 
centered at ik, where tym is small enough such that Bk does not contain any 
other ik' with k' ^ k. Take Br, with large R, such that Br contains all points 
ik- The compact set K. is given by K, = Br \ V^-i Bk and the open sets Uk are 
given by Bk \ ik, for 1 < k < N, and Uo is given by E 3 \ Br. 

Another example is a torus T 3 minus a point zo- Take a small ball B centered 
at if). Then the manifold is asymptotic euclidean with JC = T \ B and only 
one end U = B \ Iq. This is an example of an Euclidean manifold with one 
asymptotic end but non-trivial JC. More generally, given any compact manifold, 
if we subtract a finite number of points we get an asymptotically Euclidean 
manifold with multiple ends. Note that the topology of the compact core K, can 
be very complicated. 

Let us consider now Riemannian metrics on these asymptotic euclidean man- 
ifolds. Consider the manifold S — M 3 \ {0} and the metric given by (|1.17[> and 
(|1.18p . The function ip is smooth on S for any value of the constant C, however 
if C < then -0 vanished at r = —2/C and hence the metric is not defined at 
those points. That is, the metric hij is smooth on S only when C > 0. We 
have seen that S has two asymptotic ends, let us check that the metric is 
asymptotically flat (i.e. it satisfies the decay conditions <\lA\t ) at both ends Uo 
and U\. On Uq, the metric in the coordinates x 1 is clearly asymptotically flat. 
But note that in this coordinates the metric is not asymptotically flat at the 
end U\ (which, in these coordinates is represented by a neighborhood of r = 0), 
in fact the components of the metric are singular at r = 0. However, using 
a coordinate an inversion of coordinates like (jl.38p is straightforward to prove 
that the metric is asymptotically flat also at r = 0. More precisely, consider the 
coordinate transformation 




(1.39) 



In terms of this coordinates the metric has the for: 



^ = (i+£) 



(1.40) 
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We have chosen the constant factor in the coordinate transformation (|1.40[) in 
such a way that the transformation it is in fact the well known isometry of 
this metric, this choice is however not essential. The metric (jl.40|) is clearly 
asymptotically flat at the XJ\. Note that we have two energies, one for each end, 
the two are equal and given by the constant C. In this example the positivity 
of the mass is enforced purely by the global requirement of completeness of the 
metric (the energy condition is satisfied for arbitrary C). It is this condition 
that fails when C < 0. In that case the metric is defined on a manifold with 
boundary S = R 3 \B_ 2 /C: an d the metric vanished at the boundary dB_ 2 /c- In 
particular, the 2-surface dB_ 2 /c has zero area. This motivated the concept of 
"zero area singularities" introduced in ;7] , where interesting results are presented 
concerning negative energy defined on this class of singular metrics. 

In the previous example the energies at the different ends are equal. It is 
straightforward to construct an example for which the two energies are different. 
Consider the following function 

= 1 + ^+3, (1-41) 

where g is a smooth function on R 3 such that g = 0(r~ 2 ) as r — > oo and 
g(0) = a. Making the same calculation we get that the energy at one end is 
Eq = C (here we use the decay conditions on g, otherwise the function g will 
contribute to the energy at that end). But at the other end the components of 
the metric in the coordinates y l are given by 



C(l + g) 



4 



^ = { 1+ 2^) Sij > (L42) 

and hence we have that 

E 1 = C(l + a). (1.43) 

Note that in order to satisfy the energy condition (|1.20p g (and hence a) can 
not be arbitrary, we must impose the following condition on g 

Ag < 0. (1.44) 

Using (jl.44l) . the decay assumption on g and the maximum principle for the 
Laplacian (see, for example, the version of the maximum principle in the ap- 
pendix of [TB]) it is easy to prove that g > and then a > 0. 

Consider the manifold S = R 3 \ {ii}, {i 2 } with three asymptotic ends. And 
consider the function given by (this nice example was constructed in [TU]) 

* = 1 + £ + (1-45) 

2n 2r 2 

where r% and r 2 are the euclidean radius centered at the points i\ and i 2 respec- 
tively, and C\ and C 2 are constant. Note that = and hence the metric 
defined by (|1.17[) has R = 0. As before, only when Ci,C 2 > the metric is 
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smooth on S. Also, using a similar calculation as in the case of two ends it is 
not difficult to check that the metric is asymptotically flat on the three ends. 
Moreover, the energies of the different ends are given by 

Eo=Cx+C 2 , £i = d + ^^, E 2 =C 2 + ^^, (1.46) 

where L be the euclidean distance between i\ and i 2 . We see that they are 
all positive and, in general, different. These initial conditions model a head on 
collision of two black holes and they have been extensively used in numerical 
simulations of black hole collisions (see for example pQ and reference therein). 

We analyze the case of the wormhole (see which is an example of a JC 

with more complicated topology. Consider the metric on the compact manifold 
S = S 1 xS 2 given by 

7 = rf/i 2 + {dO 2 + sin 2 9dif 2 ), (1.47) 

where the coordinates ranges are — n < \i < n, and the sections /i = const are 
2-spheres. Let hij be given by 

hj = tpSij (1-48) 

where the function ij; is 

n— oo 

ip= ^ [cosh{n + 2mr}~ 1/2 . (1.49) 

n— — oo 

This function blows up at /j — 0. Hence, the metric hij is defined on S minus 
the point fj, = 0. We have seen that this is an asymptotic euclidean manifold 
with one asymptotic end. It can be proved that the metric is asymptotically 
flat at that end (see [29] for details). Also, the function ip is chosen in such a 
way that the scalar of hij curvature vanished. Moreover, the energy is given by 

oo 

E = 4^(sinh(n7r))- 1 , (1.50) 
n=l 

which is positive. 

Finally, consider the initial data for the Reissner-Nordstrom black hole given 
by a metric of the form (|1.17p with ip given by 

1 



^ =^V(<l + 2r + C)(-q + 2r + C), (1.51) 

where C and q are constant. The scalar curvature of this metric is given by 

2« 2 

Which is non-negative for any value of the constants. When C > \q\, then the 
metric is asymptotically flat with two ends Uq and U\ as in the example (|1.18|) . 
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The energy on both ends is given by E — C . The positive energy theorem applies 
to this case. If C < \q\ then the metric is singular, there is only one end Uo and 
the energy on that end is given by C . Note that in this case it is still possible 
to have positive energy < C < \q\, but the positive energy theorem does not 
apply because is a singular metric. The borderline case C = \q\ represent the 
extreme black hole. The manifold is M 3 minus a point and the metric is smooth 
on that manifold. However the metric is asymptotically fiat only at the end 
Uo, on the other end is asymptotically cylindrical. And hence this version of 
the positive energy theorem does not apply for these data. The asymptotically 
cylindrical end is a feature of all extreme black holes. For discussions on this 
kind of geometry see |17j and reference therein. 

So far, we have discussed complete manifolds without boundaries or man- 
ifolds with boundaries in which the metric is singular at the boundaries. We 
analyze now the important case of black hole boundaries. 

Black hole boundaries are defined in terms of marginally trapped surfaces. 
A marginally trapped surface is a closed 2-surface such that the outgoing null 
expansion + vanishes (more details on this important concept can be seen in 
[40]). If such surface is embedded on a space-like 3-dimensional surface, then 
the expansion + can be written in terms of the initial conditions as follows 

6+ = H - K ljS l s 3 + K, (1.53) 

where 

H = D lS \ (1.54) 

is the mean curvature of the surface. Here s l is the unit normal vector to the 
surface. For time symmetric initial data, condition + = reduces to 

H = 0. (1.55) 

Surfaces that satisfies condition (|1.55j) are called minimal surfaces, because 
(|1.55p is satisfied if and only if the first variation of the area of the surface 
vanishes. These kind of surfaces have been extensively studied in Riemannian 
geometry (see the book [31] for an introduction to the subject). We have seen 
that a marginally trapped surface on a time symmetric initial data is a minimal 
surface. That is, black hole boundaries translate, for these kind of data, into a 
pure Riemannian boundary condition. Then, we have the following corollary of 
theorem 11.1.21 

Corollary 1.2.4 (Black holes in Riemannian geometry). Let (S, hij) be a com- 
plete, asymptotically flat, Riemannian manifold with compact boundary. As- 
sume that the scalar curvature is non-negative (i.e. condition (|1.20j) ) and that 
the boundary is a minimal surface (i.e. it satisfies il.55\) . Then the energy is 
non-negative and it is zero at one end if an only if the metric is flat. 

Let as give a very simple example that illustrate this theorem. Consider 
the function rp given by (jl.181) . It is well known that the surface r — C/2 
is a minimal surface (it represents the intersection of the Schwarzschild black 
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hole event horizon with the spacelike surface t = constant in Schwarzschild 
coordinates). To verify, that we compute H for the 2-surfaces r = constant for 
the metric (|1.17p . The unit normal vector is given by 



Then we have 



Or 



Then, condition f|l .55[) is equivalent to 

°=^ + l=4-£' (L58) 

and hence for r — C/2 we have a minimal surface. Note that C must be positive 
in order to have a minimal surface. Previously we have discussed this example 
in the complete manifold, without boundaries, M 3 \ {0}. In that case corollary 
11.2.31 applies. We can also consider the same metric but in the manifold with 
boundary Rt \ Bq/2- Since we have seen that OBq/2 is a minimal surface, then 
corollary 1 1 . 2 .41 applies to that case. To emphasize the scope of this corollary, we 
slightly extend this example in the following form. Consider ip given by 

tf = (l + ^Mr), (1-59) 

where xi r ) is a function such that is x = 1 for r > C/2 and arbitrary for r < 
C/2. Corollary II .2.41 applies to this case since again the boundary is a minimal 
surface. Note that inside the minimal surface the function x is arbitrary, in 
particular it can blows up and it does not need to satisfies the energy condition. 
The corollary 11.2.31 certainly does not apply to this case. 



1.3 Linear momentum 

The total mass M defined by (|1.7[) in terms of the energy and linear momentum 
(|1.5|) - (|1.6|) represents the total amount of energy of the space-time. The first 
basic question we need to address is in what sense M is independent of the 
choice of initial conditions that describe the same space-time. That is, given a 
fixed space-time we can take different space- like surfaces on it, on each surface 
we can calculate the initial data set and hence we have a corresponding M, do 
we get the same result? We will see that the answer of that question strongly 
depend on the fall off conditions (|1.4p . 

To illustrate that, let us consider the Schwarzschild space-time. We recall 
that in the following examples the space-time is fixed and we only chose differ- 
ent space-like surfaces on it. The space-time metric is given in Schwarzschild 
coordinates (t,r s ,9,(j)) by 

ds 2 = - (l - ?p\ dt 2 + (l - 1 dr 2 + r 2 (dd 2 + sin 2 0d0 2 ). (1.60) 
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These coordinates are singular at r s = 2C and hence they do not reveal the 
global structure of the surfaces t — constant. The most direct way to see 
that these surfaces are complete 3-dimensional manifolds is using the isotropical 
radius r defined by 

l + §)\ (1.61) 
In isotropic coordinates the line element is given by 

ds 2 = -(^^j dt 2 +(^l + ^\dr 2 + d6 2 + r 2 sin 2 M/> 2 ). (1.62) 
The initial data on the slice t = constant are given by 

hij= + hi, Kij=0. (1.63) 

These are the time symmetric initial data studied in section 11.21 The linear 
momentum of these data is obviously zero, then the total mass M is equal to 
the energy E calculated in the previous section and we obtain the expected 
result M = C. 

We take another foliation of space-like surfaces. We write the metric (|1.62[) 
in the Gullstrand - Painleve coordinates (t gp ,r s ,0,(f>) (see [5S] and reference 
therein). We obtain 

ds 2 = -(1 )dt 2 gp + 2 J — dt gp dr s + dr 2 + r 2 d6 2 + r 2 sin 2 9d<j> 2 . (1.64) 

The slices t gp — constant in these coordinates have the following initial data 

/ 2m ( \ 3 
r 



hij=5ij, Kij = Uij - -SiSj ) , (1.65) 



where s l is the radial unit normal vector with respect to the fiat metric #y . We 
see that the intrinsic metric is fiat and hence the energy E is clearly zero. The 
linear momentum is also zero, because if we calculate the integral (|1.6p at an 
sphere of finite radius (note that the limit is in danger to diverge because the 
radial dependence of in (11.6510 the angular variables integrate to zero. And 
hence we obtain that for these surfaces the total mass M is zero. What happens 
is that the second fundamental form (ll.65[) does not satisfy the decay condition 
(|1.4I) since it falls off like 0(r~ 3 / 2 ). It can be proved that any initial conditions 
that satisfy (ll.4[) in the same space-time give the same total mass M. 

We consider another foliation which reveal the Lorentz transformation prop- 
erties of (E,P l ). Let (x,y, z) be the associated Cartesian coordinates of the 
isotropical coordinates (r,6,<fr), that is 



x = 7' cos sin 0, y — r sin</>cos#, z = rcos9. 



(1.66) 
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We consider the line element (11.62)) written in terms of the coordinates (t, x, y, z) 
and we perform the following change of coordinates which represents a boost in 
the z direction 

i=j- 1 (t-vz), (1.67) 
z = j- 1 (-vt + z), (1.68) 
x = x, (1.69) 
V = V, (1-70) 



where v is a constant and 7 = \/l — v 2 . Consider a surface t — constant in 
these coordinates. The intrinsic metric is given by 

h = V 4 (d.x 2 + df) + 7 ~ 2 (-iVV + tp 4 )dz 2 , (1.71) 

where 

N =Tri- (L72) 

The radius r can be written in terms of the hat coordinates as follows 



r = \J x 1 + y 2 + z 2 = \J x 2 + y 2 + 7 2 (vi+z) 2 . (1-73) 

One can check that the metric h given by (|1.71[) is asymptotically flat in the 
coordinates {x,y,z). Then, we can compute the energy of this metric and we 
obtain 

E^j^C. (1.74) 

To obtain the linear momentum we need to compute the second fundamental 
form of the slice. The calculations are long (see, for example, [20] for details), 
the final result is the following 

P x = 0, P y =0, P s = vCy- 1 . (1.75) 

Using (fl~74"l) and (fi~75|) we obtain 



M =^jE 2 - I'T'd,, =C. (1.76) 

That is, the quantities E, P l transform like a 4- vector under asymptotic Lorentz 
transformations of coordinates. 



1.4 Proof 



In section 11.21 we have presented two proofs of the positive energy theorem for 
two particular cases, for other proofs that applies to other relevant particu- 
lar cases (like spherical symmetry and weak field limit) see [9], [14], [13] and 
references therein. 
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The first general proof of the positive energy theorem was done by Schoen 
and Yau [37]. Shortly after it was followed by a proof by Witten [41] using 
completely different methods. The proof of the Penrose inequality done by 
Huisken and Illmanen in |25] (we briefly discuss this work in the following section 
II. 5p also provide a new proof of the positive energy theorem (which is based on 
an idea of Geroch [2"T] ) 

The simpler of all these proofs is, by far, Witten's one. Also it resembles 
other positivity proofs is physics: the total energy is written as a positive definite 
integral in the space. In this section we review this proof. The aim is to present 
all the relevant steps in the most elementary way. 

This proof uses, in an essential way, spinors. We refer the reader to the 
chapter in this book by Robert Geroch for an introduction to this subject. We 
will follow the notation of that chapter in this section. 

There exists various reformulations of the original proof by Witten, in this 
section we essentially follow references [34] [39] [24] [36] . 

The proof uses only spinors defined on the spacelike surface, however it is 
more transparent to begin with spinor fields in the spacetime and then, at the 
very end, to restrict them to the spacelike surface. Also, this way of constructing 
the proof easily generalize to the proof of the positivity of the energy at null 
infinity (Bondi mass) (see [36 ). 

Let [M, g a b) be a four dimensional Lorenzian manifold with connection V a . 

In this section we use the signature (H ) to be consistent with the literature 

on spinors. Unfortunately this signature gives a negative sign to the Riemannian 
metrics on spacelike surfaces used in the previous sections. 

Let Xa be an spinor field in the spacetime, the spin connection is denoted 
by Vaa', and we use the standard notation a — AA' to identify spinor indices 
with tensorial indices. 

The proof of the positive energy theorem is based on the remarkable prop- 
erties of a 2-form fl called the Nester- Witten form [41] [30], defined as follows. 
The computations of this section involve integration on different kind of surfaces 
and hence it is convenient to use differential forms instead of ordinary tensors. 
We will denote them with boldface and no indices (for an introduction to forms 
see, for example, appendix B in |40] , we will follow the notation and convention 
of this reference). 

Consider the following complex tensor 



«As/V AA'^B- 



(1.77) 



From this tensor we construct the complex 2-form f2 by 



(1.78) 



Explicitly we have 




(1.79) 



The forms used in the following are always tensor fields (usually complex) but 
they are constructed out of spinors, as in the case of $1. In order to define 
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these forms we need to antisymmetrize tensorial indices, to avoid a complicated 
notation we will always define first the tensor field in terms of spinors (as in 
equation (|1.77p ) and then define the form antisymmetrizing the tensor indices 
(as in (fTTgjt ). When there are more than two tensorial indices the explicit 
expression of the differential form (like ()1.79[) ) can be lengthy and it is not 
usually needed. The spinor has an associated (future directed) null vector 
£ a given by 

£ a = \ A \ A '. (1.80) 

Note that the $7 can not be written in terms of derivatives of pure tensors fields 
like £ a and V a . 

The strategy of the proof is the following. Consider the exterior derivative 
dtt (which is a 3-form) and integrate it on a spacelike, asymptotically flat, 3- 
surface S. Using Stoke's theorem we obtain 

^ lim & fl= I dn. (1.81) 



k 



S 



We are assuming that S is an asymptotically euclidean manifold with k asymp- 
totic ends Uk- The 2-form $7 has two important properties. The first one is 
that the left hand side of (| 1 . 81 1) gives is the total energy-momentum of a pre- 
scribed asymptotic end. The second is that the integrand of the right hand side 
is non-negative. Both properties depend on the way in which the spinor field 
A is prescribed. 

We begin with the first property. Note that the integrand in the left hand 
side of (|1.81[) is complex. But the imaginary part of fi is given by 

n-n = »V[„6] =id£, (1.82) 

where, to be consistent with our notation, we write £ for the the 1-form £ a . 
That is, the imaginary part is the exterior derivative of a 1-form and hence its 
integral over a closed 2-surface is zero. Hence the boundary integral is always 
real, for arbitrary spinors \ A . 

To prove the desired property, we need to impose fall off conditions on the 
spinor X A . Fix one arbitrary end k (from now on we will always work on that 
end, and hence we suppress the label k). Let X A be an arbitrary constant spinor, 
we require that the spinor X A satisfies on that end 

\ A = \ A + 1 A 1 7 A = 0(,-- 1 ). (1.83) 

We also assume that the partial derivatives of j A are 0(r~ 2 ) and we require 
that X A decays to zero at every other end. 

The idea is to prove that at the chosen end we have 

Pat = ^- Hm / n, (1.84) 

87T r-Hx J Sr 
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where P a = (E, Pi), with E and Pi defined by (|1.5[) - (|1.6[) . and £ a is the constant 
null vector determined by the constant spinor X by 

£ a = X A X A '. (1.85) 

Note that the boundary integral in the right hand side of (|1.84|) determines both 
the energy and the linear momentum of the end. 

To prove (|1.84[) the most important step is to prove that the value of the 
integral depends only on the constant spinor X A and not on j A . We empha- 
size, as we will see, that a naive counting of the fall behaviour of the different 
terms in $7, under the assumption (|1.83l) . does not prove this result. Using the 
decomposition (|1.83[) we write n as 

n = n + r, (1.86) 

where 
and 

r Q & = -i {Xb^aa'Ib + 7b'^aa'Xb +7s'V j4 .4'7s) , r = r [o6 ]. (1.88) 

That is, fl depends only on X A . 

We would like to prove that T = 0(r~ 3 ) and hence it does not contribute 
to the integral at infinity (|1.84[) . Consider the third term in (|1.88|) . The covari- 
ant derivative V ' aa'Ib has two terms, the first one contains partial derivatives 
of 7s which, by assumption, are Oir~ 2 ). The second term contains products 
of 7s and the connections coefficients of the space-time metric g a b evaluated 
at the asymptotic end of the spacelike surface S. These coefficients are first 
derivatives of g a b, they can be written as first derivatives of the intrinsic Rie- 
mannian metric and the second fundamental form of the surfaces and hence, 
by assumption (recall that S is asymptotically flat and hence we have the fall 
off conditions (jl.4|l ) they are 0(r~ 2 ). We conclude that "Vaa'Jb = 0(r~ 2 ) 
and hence tb'V aa'Ib = 0(r~ 3 ). We proceed in a similar way for the second 
term: since X A is constant the covariant derivative V aa' A b contains connection 
coefficients times constants and hence we have V aa'^b = 0(r~ 2 ), and then 
Tb'Vaa'Ab = 0(r~ 3 ). But using the same argument we obtain that the first 
term in (|1.88[) is 0(r -2 ) and then it can contribute to the integral. But we can 
re-write r a ;, as follows 

r a & = -i (Vbb'(7a^') -7aVbb'Aa + Jb'VaA' -Vb +7b'Va^'7s) ■ (1-89) 

The first term in (ll.89|) . which is the problematic one, contribute to T with the 
derivative of a 1-form, and hence it integrate to zero over a closed 2-surface. 
The new second term in (|1.89l) is clearly 0(r~ 3 ). We have proved that 



lim ® 11= lim d> ft. 



(1.90) 
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Note that ft is 0(r~ 2 ) and hence the integral converges. Also, the asymptotic 
value of Ct at infinity contain a combination of first derivative of the intrinsic 
metric and the second fundamental form of the surface S multiplied by the 
constants X A . It can be proved, essentially by an explicit calculation, that this 
combination is precisely P a £ a (see [41], [30] and also [4]). 

We turn to the second property of fl. Recall that the exterior derivative of 
a p-form is given by 

dn = ( P + i)v [o n 6l ... 6p] . (i.9i) 

We have 

dfl = a + f3, (1.92) 
where a and (3 are the following 3-forms 

a a bc — -«Ac'V a VbAc, a = a[ a6c ], (1.93) 

and 

Pabc = -iVaXc-VbXc P=P[abc]- (1-94) 

That is, a has second derivatives of the spinor Xa and (3 has squares of first 
derivatives of Xa- 

We compute first a. Observe that there is a commutator of covariant deriva- 
tives and hence we can replace it by the curvature tensor. However, what is 
surprising is that precisely the Einstein tensor appears. To see this, is easier to 
work with the dual of ct defined by 

*a = ^e abcd a abc . (1.95) 
We use the conmutator relations 

2V[ a V b ]Ac = —£a'B'X A bc E — £ab®a'b'c E Xe, (1.96) 

where Xabcd and $,4' b'CD are the curvature spinors. These spinors are defined 
in terms of the Riemann tensor Rated — Raa'bb'ccdd" by 

Xabcd = -^Rax'b X cy'd Y , ®abcd' — ~^Rax'b X yc ,Y d>- (1-97) 
See 33 for further details on the curvature spinors. The Einstein tensor is given 

by ' 

G ab = -6A 5ab - $ a6 , (1.98) 



where A is given by 



K=\x AB AB . (1.99) 
6 



We also use the identities 
and 



Xabc 3 = 3Ae AC , (1.100) 

Cabcd = i {f-AC^BD^A'D'f-B'C — ^AD^BC^A'C'^B'D') ■ (1.101) 
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And then we obtain 
and hence 

a = -^!^ Ge/e /^c- (1-103) 

The expressions (|1.102p and (| 1 .103(1 are pure tensorial expressions. 

To compute (3 we proceed in a similar form. We work hrst with the dual 

*P = ^e abcd (3 bcd . (1.104) 

It is important (we see later why) to split the covariant derivative V a into its 
temporal and spatial component. Let t a denote the unit timelike normal to the 
surface S and h ab is the intrinsic metric of the surface. We define the spatial 
V a derivative as 

V a = h a b V b . (1.105) 

Note that V a is not the covariant derivative D of the intrinsic metric h used in 
the previous sections, they are related by the equation 

T>AB^C — Dab^C + -^TTABC D ^D, (1.106) 

where ttabcd — ^(ab)(cd) is the spinor representation of the second funda- 
mental fundamental form of the surface. 
From equation (|1.105[) we obtain 

V Q = V a - t a t b V b . (1.107) 

We replace the derivative V a by (|1.107|) in the definition of (3 given by (|1.94[) 
and we compute the dual defined by (|1.104[) to obtain 

*P = -i^e abcd (V b \ c 'v d \ c ') + W a , (1.108) 

where 

W a = i^e abcd (t b tfV f \ c 'v d \ c + t<V V } \ c V b \ c ') . (1.109) 
Note that W a satisfies 

t a W a = 0. (1.110) 

Using the identity ([l.lOip we further decompose the first term in the right hand 
side of (ITTTU51) 

-ie abcd V b \ c 'v d \ c = V b \ B 'v BA ,\ A - V b \ A ,V AB ,X B . (1.111) 

= V C 'b\ C 'v b A ,\ A + V CB ,\ C V B ' A \ A , - V b \ A V b \ A ,, 

(1.112) 
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where in the second line we have used the spinorial identity 

tAB^CD + tBCtAD + tCACBD = 0- (1.113) 

Combining (|1.108[) and (|l.lllj) we finally obtain 

*/3 = V c , B X c 'v B A ,X A + V CB ,\ C V B ' A \ A , - V h X A V b X A , + W a . (1.114) 

We are in position now to perform the integral over S of dfl. Using (|1.92l) . 
(|1.103j) and (|1.114p we obtain 

J dfl = J (4nT ab £ b + V c >b\ C 'v b a ,\ a + V C B-\ C V B ' A \ A , - V b X A V b X A ^ t a dv, 

(1-115) 

where we have used Einstein equations 

Gab = 87TT Qb , (1.116) 

to replace the Einstein tensor by the energy-momentum tensor in the expression 
(|1.103[) for a. Note that the term W a in (| 1 . 1 14[) does not appear in the integral 
because it is orthogonal to t a (c.f. equation (|1. 11010 . 

Assume that the spinor X A has the fall-off behaviour (|1.83[) . then the identity 
(fi~84]) holds, using the Stoke's theorem ([i~8T|) (note that by (fl~83j) all the other 
boundary integrals vanish) we finally obtain the famous Witten identity 

p at = (^Tab? + V C > B \ C 'V B A ,\ A +V CB ,\ C V B ' A \ A , - V b \ A V b \ A >j t a dv, 

(1.117) 

If we assume that the energy-momentum tensor T ab satisfies the dominant en- 
ergy condition then we have 

TabCt b > 0, (1.118) 

and hence the first term in the integrand of (|1.117[) is non- negative. The last 
term in (|1.117p is also non-negative since it involves the contraction with the 
Riemannian metric (which is negative definite) and the timelike vector t AA . To 
handle the second and third term we impose on X A the following equation which 
is called the Sen- Witten equation [35] [4Tj 

V AB X A = 0. (1.119) 

Let us assume for the moment that there is a solution of this equation with the 
fall-off behaviour (ll.83[) . Then, from (I1.117[) we obtain 

Pat > 0. (1.120) 

But the constant null vector £ a is arbitrary, hence it follows that P a should 
satisfy (|1.8|) . To prove the rigidity part of theorem II. 1.11 the key ingredient is 
that E = implies, again by the identity (|1.117[) . that the spinor satisfies the 
equation 

V AB X C =Q, (1.121) 
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that is, it is covariant constant in the whole manifold. From this equation it can 
be deduced that the initial data on the surface correspond to the Minkowski 
space-time (see [32] for the details of this argument). 

It remains to discuss the solutions of equation (|1.119|) . The existence of 
solution of this equations under the required fall-off conditions (|1.83p has been 
proved in [35] [3J5] [32] ■ The main point is that equation (|1.119[) constitute an 
elliptic system of first order for the two complex components of the spinor (this 
can be easily seen using the standard definition of ellipticity for systems, see, for 
example, [16] where this specific example is discussed). And hence this equation 
can, essentially, be handled as a Poisson equation. Solutions under weak decay 
conditions on the data of equation (|1.119j) has been proved in [4] . 

Finally, let us discuss the proof of theorem 11.1.21 This was done in [3j5] [35] . 
Remarkably, the proof is very similar, the only extra ingredient is that in the 
Stoke's theorem we need to include an extra internal boundary term. This term 
has the form (see [36] ) 




In this equation + is the null expansion defined previously by (ll.53[) . The 
coefficient p' represent the ingoing null expansion on the surface, it is not im- 
portant for our purposes. The functions A and A 1 are the component of X A 
in an appropriated spinorial diad adapted to the 2-surface dB. Finally, 3 is 
a tangential differential operator to the 2-surface. It can be shown that the 
appropriate inner Dirichlet boundary for equation f)l . 1 19[) is to prescribe one 
of the component A or A 1 (but not both) (this is a consequence of the elliptic 
character of this equation, see, for example [16j for an elementary treatment of 
this). If we prescribe A 1 = on dB and use that, by hypothesis this surface sat- 
isfies + = 0, then the boundary term (|1.122[) vanished and we can proceed in 
the same way as above to prove the positivity of the energy. Note that without 
the condition + = it is not possible to make the boundary term zero. 

1.5 Further results and open problems 

In this article we have discussed only the positive energy theorem in 3 space 
dimensions. The spinorial proof presented in section 11.41 works in any dimen- 
sions (see |32] ) , however in higher dimensions the existence of an spin structure 
involves restrictions on the topology of the manifold S. In section H~4l we have 
used Weyl spinors which are well adapted to 4 spacetime dimensions. For higher 
dimensions Dirac spinors are usually used. The other proofs currently available 
|37] and |25j do not work in arbitrary high dimensions. To prove the positive 
energy theorem in all dimensions is one of the relevant open problem in this 
area. 

The positive energy theorem can be refined to incorporate other physically 
relevant parameters. For example, using a similar argument as in Witten's proof 





(1.122) 
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it is possible to prove [22] [23 that the total mass M satisfies 

M>\q\, (1.123) 

where q is the electric charge and the non-electromagnetic part of the energy 
momentum tensor must satisfies appropriated conditions. 

Recently, for axially symmetric black holes the following inequality 

M > y/\J\, (1.124) 

has been proved. Here J is the angular momentum of the black hole (see the 
review article [17] and reference therein). The equality in (|1.124|) is achieved 
only for the extreme Kerr black hole. This inequality is proved for one black 
hole, a relevant open problem is to prove it for multiple black holes. 

Another important extension is the Penrose inequality for black holes. The 
Riemannian black hole positivity theorem 11.2.41 can be generalized to include 
the area of the minimal surface, namely 

M > /A, (1.125) 

with equality only for the Schwarzschild black hole. This result was proved by 
[2"5] and [5]- The general case remains open, see the review article 27J. 

Finally, we have discussed the concept of total energy and linear momentum 
of an isolated system. It would be very desirable to have a quantity that measure 
the energy of a finite region of the spacetime. These kind of quantities are called 
quasi-local mass. For a comprehensive review on this important open problem 
see [39] ■ The following related, pure quasi- local, inequality for axially symmetric 
black holes has been recently proved 

A>8n\J\, (1.126) 

where A is the area and J is the quasi-local angular momentum of the black 
hole (see the review article [17] and reference therein). The equality in (|1.126|) 
is achieved if and only if the local geometry of the black hole is equal to the 
extreme Kerr black hole local geometry. For non-axially symmetric black holes 
it is difficult to define the quasi-local angular momentum J (sec 39]). An 
important open problem is to generalize the inequality (| 1 . 1 26|) for non-axially 
symmetric black holes (or to find suitable counter examples). 
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